Base Map erstellen
# Create a map for all radien
## Basemap mit Koordinaten von Zürich
# Zurich latitude longitude: 47.36667, 8.55
map = folium.Map(location=[47.36667, 8.55], zoom_start = 13, tiles = "cartodbpositron")In dieser Aufgabe sollen die Studierenden lernen wie man Geodaten mithilfe der Packages geopandas und folium darstellen kann. Dabei werden originale Zeckenstiche und durch eine Monte Carlo simulierte Zeckenstiche dargestellt. Informationen zu den verwendeten Datensätzen finden Sie in Table 1.
Zuerst werden die Packages geladen und die folgende Datensätze eingelesen.
| Name | Datentyp | Herkunft | Koordinatensystem | Beschreibung |
|---|---|---|---|---|
| zeckenstiche_original | gpkg | TickBite App | LV95 | Diese Zeckenstiche wurden von 2015 bis 2020 im Umkreis von Zürich über die App TickBite erfasst. Aus dem original Datensatz wurden random 100 Zeckenstiche ausgewählt. |
| zeckenstiche_simuliert | gpkg | FG Geoinformatik | LV95 | Simulierte Zeckenstiche aus Monte Carlo Simulation. Dabei wurden pro Zeckenstich 50 Runs mit 3 Distanzen (10m, 100m, 1000m) durchgeführt. |
| wald | gpkg | Swisstopo | LV95 | Dieser Datensatz beinhaltet die Waldflächen und nicht-Waldflächen in und um Zürich. |
Mithilfe von Folium wird eine Base Map erstellt, welche als Mittelpunkt die Stadt Zürich hat. Als es nächstes werden die Layer wald, zeckenstiche_original und zeckenstiche_simuliert hinzugeführt.
# Create a map for all radien
## Basemap mit Koordinaten von Zürich
# Zurich latitude longitude: 47.36667, 8.55
map = folium.Map(location=[47.36667, 8.55], zoom_start = 13, tiles = "cartodbpositron")# Wald zur Karte hinzufügen
is_wald = wald[wald["Wald_text"] == "ja"]
is_wald.to_crs(4326).explore(
m = map, # pass the map object
color = "green", # use green color on all points
tooltip = False, # hide tooltip
highlight = False, # disable highlighting
name = "Wald" # name of the layer in the map
)# Orginal Zeckenstiche zur Karte hinzufügen
zeckenstiche_original_gpd.to_crs(4326).explore(
m = map, # pass the map object
color = "red", # use red color on all points
marker_kwds = dict(radius = 1, fill = True), # make marker radius 1px with fill
tooltip = "ID", # show "ID" column in the tooltip
name = "Original Zeckenstiche" # name of the layer in the map
)# Visualisieren der Simulierten Zeckenstiche
zeckenstiche_sim_gpd.to_crs(4326).explore(
m = map,
column = "Radius",
cmap = "viridis",
legend = True,
legend_kwds = dict(colorbar = False),
marker_kwds = dict(radius = 1, fill = True),
tooltip = ["ID", "Radius", "Run"],
name = "Simulation Zeckenstiche",
)folium.LayerControl().add_to(map)
map